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Abstract: The exact partition function in ABJM theory on three-sphere can be regarded 
as a canonical partition function of a non-interacting Fermi-gas with an unconventional 
Hamiltonian. All the information on the partition function is encoded in the discrete 
spectrum of this Hamiltonian. We explain how (quantum mechanical) non-perturbative 
corrections in the Fermi-gas system appear from a spectral consideration. Basic tools in 
our analysis are a Mellin-Barnes type integral representation and a spectral zeta function. 
From a consistency with known results, we conjecture that the spectral zeta function in the 
ABJM Fermi-gas has an infinite number of “non-perturbative” poles, which are invisible 
in the semi-classical expansion of the Planck constant. We observe that these poles indeed 
appear after summing up perturbative corrections. As a consequence, the perturbative 
resummation of the spectral zeta function causes non-perturbative corrections to the grand 
canonical partition function. We also present another example associated with a spectral 
problem in topological string theory. A conjectured non-perturbative free energy on the 
resolved conifold is successfully reproduced in this framework. 
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1 Introduction 

In this work, we study non-perturbative aspects of the partition function in ABJM theory 
[1], a 3d A/" = 6 superconformal Chern-Simons-matter theory with U{N)i^ x U{N)_k gauge 
group. The ABJM theory has many important features. It is a low-energy effective theory 
on N M2-branes put on a C^/Z^ singularity. As a consequence, it has a gravity dual. In 
the ’t Hooft limit: A —>■ oo with A = N/k held finite, the dual description is type IIA 
string theory on AdS/^ x CP^. In the limit: N ^ oo with k held finite, the gravity dual is 
M-theory on AdS^ x j'Lk. We thus refer to this limit as the M-theory limit here. The 
ABJM theory is a good example to probe M-theory. 

After the seminal work by Pestun [2], a localization technique was applied to many 
supersymmetric gauge theories in various dimensions. For a recent review, see [3] and 
related articles therein. Path integrals of partition function and some BPS quantities 
finally reduce to finite dimensional matrix integrals. We emphasize that these reductions 
are exact, and no information is lost. In [4-6] (see also [7] for a review), the localization 
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was used in a wide class of gauge theories on three-sphere. An advantage of the 3d theories 


is that their matrix integrals are much simpler than those in higher dimensions. Therefore 
it is tractable to understand the large N behavior at a very quantitative level. 

Our starting point here is the exact partition function of ABJM theory [4]: 



( 1 . 1 ) 


where we follow a convention in [8]. A fundamental problem is to understand the large 
N behavior of the free energy log Zabjm(A^) ^)- This still remains non-trivial after the 
localization. In the’t Hooft limit, this was done in [8, 9] in great detail. After overcoming 
several technical issues, a recursion relation to determine the all-genus free energy was found 
[8]. An important consequence of [8] is that the genus expansion of the free energy is very 
likely Borel summable. Naively, one might expect that the all-genus Borel resummation 
reproduces the exact partition function, but this is not the case. As observed in [10], the 
Borel resummation does not agree with the exact answer. The disagreement is caused by 
so-called complex instantons [11]. The physical origin of such instantons in string theory is 
nothing but D-branes [11]. In summary, the ABJM partition function is Borel summable 
with respect to string coupling gg = 2Tr jk, but receives the non-perturbative corrections of 
form The complete large N expansion is not the genus expansion in gg but rather 

its trans-series expansion. We note that there is a framework to explore non-perturbative 
effects systematically, known as resurgence theory [12]. In matrix models and related 
topological strings, the resurgence analysis started from a pioneering work [13], based on 
[14, 15] (see also [16]). Recently, it was developed in [17-20] systematically. The resurgence 
should be useful to understand the complex instanton effects in the ABJM matrix model. 

Here we are rather interested in the M-theory limit. In an earlier important work [21], 
it was shown that the ABJM matrix model (1.1) behaves as log ^abjm ~ in the large 
N limit with fixed k. This is perfectly consistent with the M-theory expectation. However, 
the computation of the higher order 1/N corrections seems to be hard in this approach. 

A resolution to overcome this difficulty is now known. In [22], it was shown that the 
ABJM partition function (1.1) is exactly equivalent to a canonical partition function of 
an ideal Fermi-gas. The Fermi-gas approach is very powerful in the M-theory limit, and 
revealed a detailed non-perturbative structure in dual M-theory. It turned out that the 
existence of two types of instantons, i.e., worldsheet instantons and membrane instantons, 
is crucially important for the non-perturbative definition of the theory. In particular, the 
worldsheet instanton correction diverges at every physical value of the coupling, and that 
divergence is precisely canceled by the similar, but with opposite sign, divergence of the 
membrane instanton correction [23]. This pole cancellation mechanism is now observed in 
a wide class of 3d Chern-Simons gauge theories [24-27] and also in topological strings [28- 
30]. It is conceptually important since it guarantees a smooth interpolation between the 
weak coupling (type HA) regime and the strong coupling (M-theory) regime. The Fermi- 
gas formulation is also useful in ABJ theory [31], studied in [32, 33] (see also [34, 35]), in 
the topological string analysis [30, 36, 37] and in 3d mirror symmetries [38]. 
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A brief sketch. In order to make a problem clearer, let us here sketch a basic story 
in the Fermi-gas approach. A bit more detailed discussion will be reviewed in the next 
section. Throughout this paper, we focus on a generating function H(k, k) of (1.1) rather 
than the partition function itself: 

OO 

EiK, k)-.= l+Y, k^Zab3m{N, k). ( 1 . 2 ) 

N=1 

In the Fermi-gas picture, this is just the grand canonical partition function with fugacity k. 
A crucial consequence of [22] is that this grand partition function is given by a Fredholm 
determinant for density operator p describing the system. As explained in [22], the density 
operator p is of trace class, and has an infinite number of discrete eigenvalues. It is well- 
known that Fredholm determinants are entire functions and have an infinite number of 
zeros. These zeros just corresponds to the eigenvalues of p. Importantly, in the Fermi-gas 
picture, the Chern-Simons level k plays the role of the Planck constant. Therefore the 

semi-classical limit corresponds to A: —>■ 0, which is the strong coupling limit in ABJM 

theory. Thanks to this remarkable correspondence, the semi-classical analysis tells us the 
information on the membrane instantons. It is easy to see that the grand potential (the 
logarithm of the Fredholm determinant) is written as 

OO / _ 

J{K,k) :=\ogE{K,k) = ^ Z{i), (1.3) 

i=i 

where Z(s) is a spectral zeta function, defined by 

OO 

Z(s):=^A^ (1.4) 

71=0 

In the following analysis, this function plays a crucial role. Here A„ (n = 0,1, 2,...) are 
the eigenvalues of p (see (2.7)). Note that the spectral zeta function Z{s) depends on k. 
Also, note that the sum (1.4) converges only for Res > 0, but it is analytically continued 
to the whole complex plane, as in the Riemann zeta function. The relation (1.3) is already 
useful to compute the partition function for small N. In fact, using this relation, the exact 
values of Zabjm{N, k) were computed for various N and k [23, 39, 40]. To obtain the large 
N result, we need to go to the large k regime. In this regime, the grand potential takes 
the form [22] 

J{p, k) = 2^^ + M + Ac{k) -|- k), K = e^, (1.5) 

where we have introduced a chemical potential p for later convenience. In the following, 
we will use both k and p, interchangeably. From this result, one can immediately deduce 
the behavior of the free energy at large N and also the Airy function behavior first 
derived in [41]. The constant part Ac{k) is a complicated function of k, but its exact form 


- 3 - 




is known [42] (see (2.19)). The last part J7np(^, A:) is exponentially suppressed in p ^ oo.^ 
As studied in [23, 43, 44], this part turned out to be 


Jnp{f^,k)= ^ fe,rn{f^,k) eyip 



( 1 . 6 ) 


where all the coefficients k) can be computed with the help of the refined topological 

string on local x P^, in principle [28].^ Physically, the correction 0{e~‘^^) corresponds 
to the membrane instanton correction, while the correction 0(e“^^/^) to the worldsheet 
instanton correction. Note that the corrections for m 7^ 0 are invisible in the semi-classical 
expansion around k = 0 due to the exponentially suppressed correction 0(e“^^A^). These 
are understood as (quantum mechanical) non-perturbative corrections in k. 

The important point is that it is far from obvious to go from the small k result (1.3) to 
the large k result (1.5) and vice versa. One needs an analytic continuation to connect these 
two regimes. In the semi-classical limit A: —)■ 0, it was done in [22]. Also, at some special 
values of A; (A: = 1,2,4), the large p expansion (1.5) drastically simplifies, and we can 
write it down in a closed form [45, 46]. In these special cases, it is possible to analytically 
continue the large k result to the small k regime, and it reproduces all the known results. 
However this has not been done for general k so far. This problem is important because it 
is equivalent to interpolate between the small N result and the large N result, including 
all the non-perturbative corrections, directly. See also [47] for the related problem but in 
different setups and approach. 

In this paper, we give a clue to resolve this problem. Our basic tool is a Mellin-Barnes 
type integral representation.^ As in hyper geometric series, the inhnite sum (1.3) admits the 
following integral representation'^; 

rc+ioo j 

J{qi, k) = - r(.)r(-s)Z(s)e^^ (1.7) 

Jc—ioo 

where c is a constant, which must be chosen in 0 < c < 1 in the current case, as explained 
just below. To go from (1.3) to (1.7), one must be careful about the integration contour 
and the pole structure of the integrand. As discussed later, it turns out that the spectral 
zeta function does not have any poles in Re s > 0. In p < 0, we can deform the contour of 
(1.7) by adding an infinite semi-circle as shown in figure 1. Then the integral can be 
evaluated by the sum of the residues over all the poles in Res > c. If 0 < c < 1, the poles 
in this regime are at s = £ (£ G ^>0)) and thus we precisely recover the sum (1.3). On 

^The subscript “np” means the non-perturbative corrections in fi. In most of this paper, the term “non- 
perturbative” means the non-perturbative corrections in k. In this terminology, the correction for m = 0 
in (1.6) is indeed perturbative. 

^Strictly speaking, the grand potential receives an additional contribution to the topological string 
prediction. This was first observed in [23] , and is now understood as a contribution from a generalization of 
Jacobi theta function [30, 45, 46]. Here we do not care about it because it is not important in the following 
analysis. The important fact is that this contribution does not change the form (1.6). 

®We thank Marcos Marino for telling us about a powerfulness of the Mellin-Barnes representation. 

“^This is essentially an inverse Laplace transform. For a relationship between spectral zeta functions and 
spectral determinants, see [48], for instance. 
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Figure 1. (Left) (Expected) pole structure of Z(s) at fc = 6 is sketched. The blue “x” show 
“perturbative” poles while the red “o” show “non-perturbative” poles. In this case, the non- 
perturbative poles appear at s = —2/3, —4/3, —2,.... The non-perturbative poles at s = —2, —4,... 
are overlapped with the perturbative poles. At these point, the pole cancellation [23] in the grand 
potential occurs. (Right) Contour deformations of the Mellin-Barnes integral (1.7). For /i < 0 
(/r > 0), we can add the semi-circle C+ (C-). We also show the pole structure of the integrand in 
(1.7) at /c = 6. To pick up the poles correctly, we have to choose c such that 0 < c < 1. 

the other hand, for fi > 0 we can deform the contour by adding the opposite semi-circle 
C- as in figure 1. In this case, one needs the information of the poles in Res < c. In 
this regime, Z{s) may have non-trivial poles, and the problem is highly non-trivial. From 
the Mellin-Barnes integral point of view, the large // expansion (1.5) should be understood 
as the sum of the residues over all the poles in Re s < c. More explicitly, the large 
expansion, from (1.7), should be given by 

jU,k) = - y Resr(s)r(-s)Z(s)e"^, Ul^oo). (1.8) 

All the poles 
in Re s<c 

As we will see in the next section, the polynomial part in (1.5) indeed comes from the residue 
at s = 0. The contributions from the poles in Res < 0 are exponentially suppressed in 
/r —>■ oo. Therefore we conclude that to reproduce the result (1.5) with (1.6), the integrand 
of (1.7) must have the poles (only) at 

A.7TI 

s = -2e-—, £,m = 0,1,2,... . (1.9) 

k 

This is a consistency requirement from the known results, and must be confirmed by a 
direct analysis. Since the gamma function r(±s) has poles only at s = 0,t1)T2, ■■■, 
respectively, we conclude that most of the poles in (1.9) must come from the spectral zeta 
function Z{s). We note that the poles for m yl 0 go to the infinity in the limit A: —>■ 0, and 
they are invisible in the semi-classical expansion. In this sense, we call these poles non- 
perturbative poles here.^ In the semi-classical expansion, only the poles with m = 0 appear. 

®A concept of non-perturbative poles is also found in a similar but different context [34, 52]. We note 
that here we consider the small k expansion while there the large k (or small Qs) expansion was considered. 
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Spectral Problem 
(Fredholm Integral Equation) 

I 

Spectral Zeta Function 

I 

Mellin-Bames Integral 
(Grand Potential) 

^ N 


Small K. Expansion 
Small N Result 


Earge Hi Expansion 
Earge N Result 


Figure 2. Basic flow of our analysis. This approach is widely applicable to other examples with a 
well-defined spectral problem. 

At present we have no way to determine the complete pole structure of Z(s), but in this 
paper we give an evidence that it has at least a non-perturbative pole at s = —A/k. The 
Mellin-Barnes representation states that this pole produces the correction which is 

nothing but the leading non-perturbative correction to the grand potential. This idea is 
widely applicable to other examples. In fact, in addition to the ABJM Fermi-gas, we also 
pick up another example, in which this mechanism works. This example is related to a 
spectral problem in topological string theory [30, 49, 50]. We show that a conjectured non- 
perturbative correction to the free energy [28] on the resolved conifold is indeed reproduced 
in this approach. 

Our basic strategy in this paper is schematically shown in figure 2. The Mellin-Barnes 
type representation (1.7) is very involving. It is clear that all the information on the 
(grand) partition function is encoded in the spectral zeta function in the whole regime of 
H (or N). In particular, to extract the large /x result, it is important to understand the 
pole structure of Z{s). We observe that the resummation of the semi-classical expansion of 
Z{s) produces the non-perturbative pole at s = —Ajk. This is conceptually interesting. We 
only consider the semi-classical analysis on the spectral zeta function. Nevertheless we can 
explain the appearance of the non-perturbative correction to the grand potential.® This is 
because that the Mellin-Barnes integration does not commute with the WKB infinite series. 
The perturbative resummation of the spectral zeta function causes the non-perturbative 
corrections to the grand potential. This is a main message of this paper. 

®There is, however, a possibility that the spectral zeta function itself receives non-perturbative correc¬ 
tions. If so, these corrections might also have the poles (1.9). We have no evidence of the existence of such 
corrections, but cannot exclude this possibility in this work. We thank Marcos Marino for pointing out this 
issue. 
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Outline. The organization of this paper is as follows. In section 2, we briefly review 
the Fermi-gas approach proposed in [22]. In particular, we see that the Mellin-Barnes 
representation (1.7) is indeed powerful to know the grand potential in the large /r limit. In 
section 3, we propose a toy model, which has many common features to the ABJM Fermi- 
gas system. Through this example, we can learn how the non-perturbative corrections to 
the grand potential appear via the Mellin-Barnes representation. In section 4, we proceed 
to the analysis of the ABJM spectral zeta function. First, we develop an efficient way 
to compute the semi-classical expansion systematically. Then, using this result, we show 
that the spectral zeta function indeed has the non-perturbative pole at s = —4//c. In 
section 5, we give another example, which is closely related to a spectral problem for the 
topological string theory [30]. We see that our prescription here indeed reproduces the 
non-perturbative free energy of the resolved conifold. In section 6, we give conclusions. 
Appendix A is a detailed computation of the Wigner transform, which is used in section 4. 
In appendix B, we give a few results on the differential operators, which are useful to 
compute the quantum corrections to the grand potential and to the spectral zeta function. 


2 Review of the ABJM Fermi-gas 


Let us quickly review the Fermi-gas formulation in [22]. We start with the exact partition 
function (1.1). The crucial step is to rewrite this matrix integral as the following form 

1 /•oo ^ 

ZABm{N,k) = ( 2 . 1 ) 

^ o J — OO -1 


where 


p(xi,X2) = {xi\p\x2) = 


o-gSn 

1 1 


i=l 


1 


1 


( 2 . 2 ) 


27:k (2 cosh f) 1/2 (2 cosh (2 cosh f )i/2 ' 

Importantly, the expression (2.1) can be regarded as the canonical partition function of 
an ideal Fermi-gas described by the unconventional density operator p. The quantum 
Hamiltonian of this system is given by 


H 


-'—1 I 

= p = 2 cosh — 


x\ 


1/2 


2 cosh ^ 
2 


2 cosh — j 


(2.3) 


where x and p are the canonical variables satisfying the commutation relation 


[x,p]=ih, h = 27rk. (2.4) 

Therefore in this formulation, the Chern-Simons level k plays the role of the Planck con¬ 
stant. The semi-classical limit is A: —>■ 0, which is a very strong coupling limit in ABJM 
theory. Since the operator p is a non-negative, Hermitian and Hilbert-Schmidt operator, it 
has positive discrete spectra Aq > Ai > A2 > • • • > 0. As already noted in the introduction, 
for the canonical partition function (2.1), the grand partition function (1.2) is given by the 
Fredholm determinant: 

OO OO 

H(k,/ c) = det(l-h k/ 5) = JJ(1«;A„) = ]^(1e'"”®"), = (2.5) 

n=0 n=0 
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where En is the n-th energy eigenvalue of the system. The eigenvalue problem is now given 
by 

p\(t>n) = K\<i>n) , (n = 0,1,2 ,...). (2.6) 

As shown in [22], this can be recast in a homogeneous Fredholm integral equation of the 
second kind in the coordinate representation: 


f 


dx' p{x,x')4)n{x') = \n4>nix) 


(2.7) 


One can easily rewrite (2.5) as the form in (1.3). Therefore if one knows all the eigenvalues 
in (2.7), one can compute the spectral zeta function, in principle. However this is not easy 
for generic k. 

The first step is to consider the semi-classical approximation, as in quantum mechanics. 
Let us consider the semi-classical limit k ^ 0. At the leading approximation, we can treat 
X and p as commuting variables. The classical Hamiltonian is thus given by 

_ ^2 cosh ^2 cosh . (2.8) 


Then the spectral zeta function for s > 0 is easily computed by the phase space integral 

[ 22 ]: 



1 r^{s/2) 
4^ r2(s) 


(2.9) 


This is easily analytically continued to complex s. Substituting this result into (1.3), one 
gets an analytic form of the grand potential in terms of the hypergeometric functions [22]. 
Then, it is possible to analytically continue the grand potential to the large p regime. Here 
we use the Mellin-Barnes integral (1.7) to extract the large p expansion directly from the 
spectral zeta function.^ The classical grand potential is represented by 




1 


L 


c+ioo 




r(s) ’ 


( 2 . 10 ) 


where as explained in the previous section, we have to choose c such that 0 < c < 1. It is 
obvious that the integrand in (2.10) has poles at s = —2n (n = 0,1,2,...) in the Re s < c 
region. For the pole at s = 0, one obtains the contribution 




s/i _ 


2p^ 


P 

+ w + 


2C(3) 


TT^ 


( 2 . 11 ) 


Similarly, from the poles at s = —2n (n > 1), the exponentially suppressed corrections 
with are obtained. We conclude that the classical grand potential has the following 

large p expansion: 




4/i 4 2 

H —^ H —^ ~ ^ 


TT^ 


TT^ 


(p^oo). (2.12) 


^This approach has already been discussed by Marcos Marino in his review article [51]. 
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This, of course, reproduces the result in [22]. We stress that to obtain this result, one does 
not have to perform the sum (1.3). The Mellin-Barnes integral (1.7) allows us to know 
the grand potential in the large /r regime directly from Z{s). As explained in section 4 in 
detail, it is possible to compute the semi-classical corrections to the spectral zeta and the 
grand potential systematically. These quantities admit the WKB expansions 

^ OO - oo 

Z^nkb{s) = (2.13) 

n=0 n=0 


As in (4.15), the first correction Z^^^s) is given by 

zC)(,) = =!£!(TAz(0)(,), 

^ ^ 96 (1 + s) ^ ^ 


(2.14) 


Using the Mellin-Barnes representation again, one easily finds 

Pushing the same computation in the higher corrections, we observe that these corrections 
start from the constant term, i.e. 


= 0{l), n>2, 


(2.16) 


Thus the WKB expansion of the grand potential is generically given by 
o 3 / 1 i. \ 

JwKB{fJ',k) = + — j + Ac{k)+ + bi{k)n + Ci{k))e (2.17) 

where the constant term has the non-trivial corrections in all orders in small k expansion. 
It is not easy to determine its exact form only from the semi-classical analysis, but the 
all-order prediction is now known [42] because this part corresponds to the constant map 
contribution in the topological strings. The all-order small k expansion is given by 

■4c(^ = ^ ^ E (2,18) 

n=l ' ' 

where Bn is the Bernoulli number. Obviously this sum is a divergent series, but Borel 
summable for A; G M. After resumming it, one finally obtains the following integral repre¬ 
sentation [24]: 

In the above semi-classical analysis, only the correction with m = 0 in (1.6) appears. 
In fact, at each order in the WKB expansion, we encounter only the poles at s = —2n 
(n = 0,1, 2,...) for Re s < c. Physically, these corrections are interpreted as the membrane 
instanton corrections [22]. We note that all of the coefficients a^^k), hi{k) and ci{k) in (2.17) 
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are exactly predicted from the refined topological string on local x P^ [28], in principle. 
The other corrections in (1.6) appear as non-perturbative corrections in k. These are due 
to the worldsheet instanton corrections and bound states of the two kinds of instantons. 
The main purpose in this work is to explore a mechanism how these non-perturbative 
corrections are produced. 

Remarkably, if introducing the following redefined chemical potential 

2i OG 

PefT := lip (2.20) 

i=i 


then the grand potential (1.5) drastically simplifies® [44] 






\3k “^247 


IJ.eS + Ac{k) + JM2{fJ-eS, k) + Jws(/^efr, k), 


( 2 . 21 ) 


where 

OO 

JynilJ^^k) = ^(6£(A:)pQ( fc))e“^^^, 
i=i 

OO 

JwsilJ;k) = ^ dm{k)e 

m=l 


( 2 . 22 ) 


In the remaining sections, we see that the resummation of the perturbative WKB sum (2.13) 
of Z{s) causes the non-perturbative poles. These poles are sources of the non-perturbative 
corrections to the grand potential in (1.6) or (2.22). 


3 A toy model 


In order to understand the essence, let us start with a very simple toy model. We consider 
the following exact spectral zeta function: 


z*°y(s) 


I fks ks\ 

l^T’Ty ’ 


(3.1) 


where B(x, y) = r(x)r(y)/r(x + y) is the Euler beta function. This spectral zeta function 
admits a semi-classical expansion® around k = 0 


Z^°y(s) 


2 _ C(3)^ 2 

n'^s'^k 48 Ibvr^ 




5120 


k^ + Oik"^). 


(3.2) 


Note that Z*°^(s) does not receive any non-perturbative corrections in k. We want to 
understand the large p behavior of the grand potential in this toy model. First, it is easy 
to see that has infinite poles at 

Am 

s = ——, m = 0,1,2,.... (3.3) 

k 

®As noted in the introduction, strictly, we need to subtract an “oscillatory part” in [23]. Here we do not 
care about it. However, it will be important when one wants to compare the coefficients at higher orders. 

^The WKB expansion (3.2) is slightly different from (2.13) for ABJM spectral zeta. The expansion (3.2) 
is rather similar to the so-called Nf matrix model in Nf 0 [24, 53, 54]. 
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There are no singularities except for these points on the whole complex plane. Using 
the Mellin-Barnes representation (1.7), one can know the large ^ behavior of the grand 
potential. In addition to the above poles, the integrand of the Mellin-Barnes integral has 
poles at s = —n (n = 1, 2,...) in the region Res < c. The pole at s = 0 determines the 
leading behavior in the large /r limit. It is easy to find 


where are subleading exponentially suppressed corrections in 

/r —)• oo. We separate these two contributions because their origins are slightly different. 
The “membrane instanton” correction J^{iJL,k) comes from the poles at s = —re. The 
terminology is just an analogy with the ABJM Fermi-gas. We find that it is generically 
given by 


JmI= Y1 

n=l 


(-ir p 

(27rre)2 



(3.5) 


On the other hand, the “worldsheet instanton” correction comes from the non- 

perturbative poles s = —Am/k. After a simple computation, one finds 





(3.6) 


The result (3.4) with (3.5) and (3.6) is the complete large // expansion of the grand potential 
for the toy model. In this example, we can completely control all the corrections. 

Let us remark some of important lessons that we should learn from this toy model. 

Pole cancellation. We notice that both {fJ-, k) and 77^g (/r, k) diverge at rational 
values of k. As an example, let us consider the limit /c —> 4. In this limit, the membrane 
instanton part is given by 


71=0 


2(-l)’" 
n^'7r'^{k — 4) 


2 n\ (- 1 )" 

re J re^vr^ 



{Hn-H2n) + 0{k-A) 


e-n/i^ 


(3.7) 


where Hn is the re-th harmonic number. The first term is divergent. Similarly, the world- 
sheet instanton part is given by 


Jwg(/x,A:) = ^ 

71=0 

The important point is that the sum of these two corrections is finite 



This is, of course, not accidental. At the level of the spectral zeta function or the Mellin- 
Barnes representation, there is no problem in the limit k ^ A. It gives finite values for any 


2(-l)’" /2re 
(/c — 4) I re 


re^TT^ 


+ 


(-ir 

re^TT^ 


2re 

re 


M + i 

2 re 


+ C>(A:-4) 


^-n/i 


(3.8) 
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k. What happens in /c —>■ 4? In this limit, the non-perturbative poles s = —injk “collide” 
with the perturbative poles s = —n. At A: = 4, both poles degenerate, and the orders of 
the poles become higher. One can easily check that the spectral zeta function at A: = 4, 

= (»,«), (3,10) 

correctly reproduces the mixed corrections (3.9). 

Non-perturbative effects. The toy model here is already involving. At the level of 
the spectral zeta function, there are no non-perturbative corrections in k. Everything 
can be understood only in the perturbative analysis. Nevertheless, the grand potential 
receives the non-perturbative correction J^^{lJ^,k) in the large /r regime. The source of 
this correction is the non-perturbative poles at s = —Am/k {m = 1,2,...). We stress 
that these poles are never visible in the semi-classical expansion around A; = 0. In fact, 
using the WKB expansion (3.2), one encounters, at each order, only the perturbative poles 
of the integrand in the Mellin-Barnes representation. The non-perturbative poles appear 
after resumming all the perturbative correetions to Z*°y(s). In the next section, we will 
give a strong evidence that the similar mechanism also works in the ABJM Fermi-gas: The 
perturbative resummation induces non-perturbative poles. 

4 Spectral zeta function in the ABJM Fermi-gas 

Let us proceed to the ABJM Fermi-gas. Unfortunately, unlike the toy model in the previous 
section, the exact spectral zeta function in this model is not known.Therefore we cannot 
apply the method in the previous section immediately. The situation is very limited. To 
explore its analytic property, we first consider the semi-classical expansion of Z{s), as in 
(2.13). For s = i, (^ = 1,2,...), there is a powerful way to compute the semi-classical 
expansion of Z(i) by using the TBA equations [43]. Here, however, we take another 
strategy. We use the Wigner transform in quantum mechanics. This method was already 
considered in the original paper on the ABJM Fermi-gas [22], but we develop a bit more 
efficient way. This approach is useful in the case where the TBA description is not known. 
We next extrapolate the semi-classical expansion to the finite-A; regime by using the Fade 
approximation, and observe the pole structure of Z{s). 

4.1 Semi-classical analysis: Wigner transform 

For a quantum mechanical operator A, the Wigner transform is defined by 

Aw{x,p):=J dx'e~{x-—\A\x + —), (4.1) 

where recall that h = 2TTk. By definition, the diagonal element of A in the coordinate 
representation is given by 

{x\A\x) = J -^Aw{x,p). (4.2) 

should note that, at fc = 1, 2, 4, the exact spectrum can be computed by solving an exact quantization 
condition [29, 45, 46, 55, 56]. It would be interesting to analyze the exact spectral zeta on the complex 
plane for these special cases. 
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Thus the trace is given by 


Tr^ = 


/ OO 

dx(x|^| 

-OO 


X) = 


dxdp 

27rh 


Aw{x,p). 


(4.3) 


We hrst apply the Wigner transform to the inverse of density operator 


O := p ^ = ( 2 cosh ^ 


1/2 


P 


2 cosh (r 2 cosh — 


X 


1/2 


(4.4) 


As derived in appendix A, the Wigner transform of O is exactly given by 

1/2 


Owix,p) = 4 cosh ^ (^cosh^ ^ — sin^ 


We want to compute 


2 V 


Z{s) = Trp* = Tre"^" = TrO 


4 / 


(4.5) 


(4.6) 


We use a general argument in [22]. For a given function /, the Wigner transform of f{0) 
is given by 




r=0 


where 


Gr ■= [{O - Ow)^]w- 


(4.7) 


(4.8) 


In computing this quantity, one encounters the Wigner transform of operator products. 
The Wigner transform of a product of two operators is computed by 


(A • .B)w = * -Bw, 


(4.9) 


where the Moyal product * is defined by 
A-k B := A{x,p) exp 

OO n 

= EE(-i) 


ih ,k- —I ^ 

— {dxdp - dpdx) 


B{x,p) 


n=0 m=0 


")4 (I) 


(4.10) 


In our purpose, we set f{x) = x In this case, the Wigner transform (4.7) becomes 

(a-*)w = f; (4.11) 

r =0 

where (s)^ is the Pochhammer symbol. Using the exact relation (4.5), one can compute the 
WKB expansion of this formula up to any order, in principle. Then after the integration 
over {x,p), one obtains Z{s). Note that the integral converges only for Res > 0. After the 
integration, one can do the analytic continuation. 

It is complicated to compute the integral of (4.11) for general s. Practically, it is 
sufficient to compute it for s € Z>o. Prom the lower order information (see (4.15)), one can 
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guess a form of higher order correction with finite unknown parameters. These parameters 
are fixed only by the result for finite s G Z>o. If the ansatz is correct, the obtained result 
must reproduce the corrections for other values of s. In this way, one can fix each coefficient 
in (4.15) order by order. Equivalently, it is convenient to use an interesting observation 
in [26]. Recall that the WKB expansion of the grand potential is given by (2.13). The 
quantum correction (fi) can be constructed by acting a non-trivial differential operator 
to the classical one: 

= n = l,2,..., (4.12) 

where is a differential operator of Its explicit form up to n = 2 is found in [26]. 
Using this observation, we have indeed fixed the differential operator up to n = 10. The 
result up to n = 4 is given in appendix B. The obtained grand potential in this way is in 
precise agreement with the one from the TBA approach [43].^^ 

Now let us consider the WKB expansion of the spectral zeta function, given by (2.13). 
As already computed in section 2, the classical part is explicitly given by (2.9). Combining 
(1.3) and (4.12), one easily obtains the relation 

^(0)(^)pW(e^M) = ^(n)(^)e^M, ^ = 1,2,.... (4.13) 


This means that the correction {£) is obtained by replacing 9^ in by £. It is very 
natural that the correction is analytically continued by 

^(0)(5)p(n)(eSM) = ^(n)(^)eSM, ^ g c. (4.14) 

Using the explicit form of in appendix B, the WKB expansion of Z{s) is given by 


■^wkb('S) = 


Z(0)(i 

k~ 


5^(1 - s) 


s^{l-s){l6 + 27s + 7s^) 




96(1 + s) 


92160(1 + s)(3 + s) 


5^(1 - s) (256 +2560s + 3216^2+ 1649 s3 + 376s^ +31s®), ' 


(4.15) 


61931520(1 + s)(3 + s)(5 + s) 

We use this expansion in order to explore the pole structure of Z(s) in the next subsection. 


4.2 Pole structure and non-perturbative effects 

In this subsection, we explore the pole structure of Z{s) for finite k. Since we do not know 
the all-order coefficients of the WKB expansion, we cannot extrapolate it to the finite-A: 
regime, directly. Here we apply the Fade approximation to the WKB expansion (4.15). An 
advantage of the Fade approximation is that it captures the pole structure of the original 
function. We see that the extrapolation of the Fade approximant of (4.15) to finite k 
strongly implies that Z(s) has a pole at s = —4:1k, which is of course invisible as long as 
one considers only each term in (4.15). 

We thank Kazumi Okuyama and Masazumi Honda for sharing their results on the higher order differ¬ 
ential operators in another project [57]. 
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Let us first see that the spectral zeta function Z{s) does not have any poles in Re s > 0. 
In this region, we can use the definition (1.4). We use an inequality 




OO 


oo 


OO 


y~]e 

n=Q 


< ^ |e 

n=0 


n=0 


(4.16) 


where x := Res and = — logA^ is the eigenvalues of the quantum Hamiltonian H. 
Since £■„ > 0 for Vn G Z>o, the rightmost equation monotonically decreases w.r.t. x. We 
conclude that the rightmost equation is finite for x > 0, and thus Z{s) has no poles in 

X > 0. 

To see the poles in Re s <0, let us consider the Fade approximation of a given function 

OO 

fix) = 

n=0 


We denote its Fade approximant with numerator order M and denominator order N by 


[M/N]f{x) 


ao + oix + • • • qmx^ 
1 + bix + • • • b^x^ 


(4.18) 


The M + N + 1 coefficients {oq, ...,om; bi,... are uniquely fixed by requiring 


/(x) - [M/N]f{x) = C)(x^+^+^). 


(4.19) 


The Fade approximant is the most general form in the ratonal-type approximation. It, of 
course, includes the original Taylor expansion. It is well-known that the Fade approximant 
captures the global structure of the original function, compared to the Taylor expansion. 


Test in the toy model. Before proceeding to the ABJM spectral zeta, let us test a 
validity of this method by using the toy model in the previous section. We start with the 
WKB expansion (3.2). For convenience, we denote this expansion by 








ytoy 

^WKB 


(s;fc) = 1 - 




TT S 


96 


+ 


32 


10240 


+ Oik^). 


(4.20) 


Then, we apply the Fade approximation to Z^^g(s; k). 


We also introduce a notation 


ytoy,[M/N] 
^Pade 


{S]k) 


n'^s^k 


[M/N]^toy {s-,k). 

^WKB 


(4.21) 


What we show here is to compare this approximant with the exact function (3.1) for finite 
k. In figure 3, we show the behavior for k = 6,8. The red solid lines are 

the approximant Z^^^f^^\s',k), while the blue dashed ones are the exact Z^°^{s',k). Both 
graphs show good agreements in the regime —Ajk < s < 0. 


- 15 - 












Figure 3. This figure shows the Fade approximant (Red solid line) and the exact 

function Z^°^{s\k) (Blue dashed line) at fc = 6 (Left) and k = 8 (Right) as functions of s. The 
Fade approximant captures a global structure of the original function (at least to the nearest pole 
from the origin). 




Figure 4. We show the Fade approximant (s; fc) at /c = 5 (Left) and k = 6 (Right) as a 

function of s. It is observed that they diverge at s « —4/fe. 


Application to ABJM. Now let us apply the Fade approximation to the spectral zeta 
function in the ABJM case. As in the toy model, we first define 


■Z’wkb(s; k) 


fcZwKB(s) 

Z(o)(s) 


1 + 


g^(l - s) 
96(1 + s) 


(tt/c)^ H-. 


(4.22) 


We use the WKB expansion up to 0{k'^^) and compute the Fade approximant [10/10]^^^^. 
We also denote 

4T'(o;^ = (4.23) 

In figure 4, we show the behavior of (s; k) for k = 5,6. These graphs strongly suggest 

that Z{s) has a pole at s = —Ajk. Using the Mellin-Barnes representation, this pole leads 

_ 4fj, 

to the non-perturbative correction of order e~~. The result is completely consistent with 
the leading worldsheet instanton correction. 


Non-perturbative corrections. It is clear that the coefficients of the non-perturbative 
corrections are related to the residues of the non-perturbative poles. If the spectral zeta 
function behaves near the non-perturbative pole at s = —Am/k {m = 1,2 ,...) as 


hm Z[s) = --—— 

s^-im/k s -b 4m/k 


+ o{i), 


(4.24) 
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then the Mellin-Barnes representation states that the grand potential receives the correction 


4m 


zir, {k) CSC 


(T) 


Amjj, 

'~ir~ 


(4.25) 


where we have assumed that this pole do not overlap with any other poles. To reproduce 
the known result^^ in [23] the first residue z\ ^\k), for example, should take the form 


z\ ^\k) = — cot 

TTK 



(4.26) 


So far, we have no principle to determine such residues analytically. It would be very 
interesting to understand the analyticity of Z{s) in more detail. 


4.3 A comment on non-perturbative correction to spectral zeta function 

In the previous subsection, we saw that the WKB perturbative resummation produces the 
non-perturbative pole that causes the leading non-perturbative correction to the grand 
potential. We must note that there is a possibility that the spectral zeta function itself 
receives non-perturbative corrections of form This type of correction is invisible in 

the semi-classical analysis. Let us give a comment on this possibility. We first note that 
at s = 1, 2 the spectral zeta function is exactly computed in [58]: 

1 I r°° f 

Z{1) = Z{2) = ^ / dt -5-. (4.27) 

dk 47r2 Jq cosh^ t sinh kt 

It is easy to see that these results do not receive the non-perturbative correction in the 
small k limit. Also, if plugging s = —2, —4 into (4.22) (and using the result in appendix B), 
we get the following expansions: 


Zwkb{s 

Z^nkb{s 


-2-k) 

-A-k) 


(TTky _ {T^kf {irkf oik^^) 

8 384 46080 10321920 ^ 

^{T^kf h{'Kkf _ {Tlkf {Tlkf 

18 216 1296 72576 ^ 


It is easy to guess that these are reproduced by 


(4.28) 


■Z’wkb('S 
■^wkb('S 

Hence for s = —2, —4, the WKB sums can be performed exactly.^^ Obviously these do 
not have the non-perturbative corrections in k. All of these results seem to imply that 

^^Recall that to compare the known result one should be careful about the contribution from the gen¬ 
eralized theta function [30, 45, 46]. At this order, however, there is no such a contribution, and one can 
compare the result with the topological string one. 

^®This is not accidental. In fact, the function Zwkb(s = —2£-,k) are related to the function ai{k) in 
(2.17). This is a consequence of the property (4.12). We also note that at s = —2£, the classical part 
is divergent but ZwKB{s',k) is finite. 


-2;k) = cos—, 

-4; fc) = - (4 -|- 5 cos 7rk). 
9 


(4.29) 
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the spectral zeta function does not receive the non-perturbative corrections. However, we 
cannot conclude it at present. We need a more detailed analysis. If the spectral zeta 
function has the non-perturbative correction, these corrections might also have the non- 
perturbative poles in (1.9), in principle. We hope to report a resolution of this issue in a 
future. 

5 Another example: topological string on the resolved conifold 

In this section, we give another example, for which our procedure works. In particular, 
we will see that the topological string free energy on the resolved conifold is precisely 
reproduced by our prescription in a non-perturbative manner. This example is closely 
related to a quantum spectral problem in the topological string theory [30, 50]. The setup 
here is explained in [30] in detail. 

5.1 General framework 

In [30], a new perspective to formulate the topological strings was proposed, based on 
an earlier work [49] and the Fermi-gas formulation [22]. The key idea is to associate an 
operator of a quantized mirror curve with a spectral problem. Our starting point is a mirror 
curve describing a toric Calabi-Yau threefold. Here we restrict ourselves to the following 
simple case: 

W(e",e^>) = O^,„(x,p) + u = 0, (5.1) 

where Om,n{x,p) is defined by 

0^,n{x,p) := e" + e^’ + . (5.2) 

As explained in [36, 37], this mirror curve describes the anti-canonical bundle of the 
weighted projective space P(l, m,n). Now we want to “quantize” this curve. Following the 
prescription in [49], the quantized mirror curve is given by 

(e^ + eP + e-™*-”P + {i)lT) =0, (5.3) 

where x and p satisfy the canonical commutation relation (2.4). The wavefunction T(x) = 
(xlT) describes a brane that probes the geometry [49]. As conjectured in [30] and proved 
in [36], the inverse operator 

Pm^n := = (e" + eP + g—(m ,n > 0), (5.4) 

is positive-definite and of trace class. This means that the operator pm,n has positive and 
discrete eigenvalues, as in the ABJM Fermi-gas. In [30], an explicit formula for a spectral 
determinant of this operator was conjectured (for some special cases). In [36], it was shown 
that this operator has a good representation in certain coordinate. This representation is 
useful to reformulate the topological strings as matrix models [37]. Here we focus on the 
spectral zeta function Zm,n{s) of the operator pm,n- 
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As in sections 2 and 4, we consider the semi-classical limit h —>• 0. As before, the 
spectral zeta function and the grand potential have the following WKB expansion: 


- LXJ ^ LXJ 

e=o ^ £=o 

At the classical level, the spectral zeta function is computed by 

’dxdp 1 _ r()r(^)r() 


Z^U^) = I - 


2-k qP + Q-mx-npy 2'n{m + n + l)V{s) 

Therefore the classical grand potential is given by 


27r 27ri 


,r(-s)r 


m + n -|- 1 


ms 


m + n + 1 


ns 


m + n -|- 1 


(5.5) 


(5.6) 




(5.7) 


In the large p limit, one finds 

, (m-|-n + l)^ q Tr(m‘^ + mn + n'^ + m + n + 1) 

4?,U/^) = ^-TTyrz- 


+ 


127rmn 12mn 

(m + l)(n + l)(m -|- n)C(3) 


.. , m-fn + l .. m-fn+1 

+ 0{e - '^,e " 




(5.8) 


2TTmn{m + n -|- 1) 

where exponentially suppressed corrections come from the residues of the poles in s < 0. 

Next we consider the Wigner transform to compute the quantum corrections. As in 
the previous section, the Wigner transform of the operator Om,n can be computed exactly, 
and it takes the very simple form 


(Om,n)w(x,p) = e" + e^'+ e-“P. 


(5.9) 


The Wigner transform is just the same form as the original function! One can confirm 
this result by a direct computation along the line in appendix A. Then using (4.11), one 
can compute the quantum corrections to Zm,n{s) systematically. As in (4.12), it is useful 
to fix the differential operator acting on the grand potential. In appendix B, we present 
the explicit forms of the first three corrections. Using these operators, one finds the WKB 
expansion of Zm,n{s)'- 


z^Ti^) = 


ZJm,.n, \ * ^ 


h 


1 - 


mns 


24(m + n + 1) 




+ 


mn{m + l)(n + l)(m -|- n)s^ 


+ 


7m?n^s'^ 


2880(m-k n-M)2 5760(m + n + l)2 

and also the grand potential in the large p limit: 

:7™(/«, ^) = + Am,n{K) + 

where • • • denotes the exponentially suppressed corrections and 

(m + n + 1)^ 


+ 0{h^ 


Bm,niK) = 


Anmnh ’ 

7r(m^ -|- mn + -|- m + n + 1) (m -|- re + 1)H 


12mnh 


487r 


(5.10) 


(5.11) 


(5.12) 
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Finally we conjecture the form of Am,n{K) with the full quantum corrections: 


Am,niK) 



A, 


(m + n + 1)^ 


vr 


(5.13) 


where Ac{k) is just the same function appearing in (2.19). Using the differential operators 
in appendix B, this conjecture is checked up to h^. It is not easy to determine exact forms 
of the exponentially suppressed corrections. All of these results are valid for any m > 0 and 
n > 0, for which the spectral problem for pm,n is well-defined. In particular, for m = n = 1, 
the corresponding geometry is local P^. In this case, all the results here reproduce the ones 
in [30]. In the next subsection, we analytically continue them to n < 0. In this case, 
the spectral problem is no longer well-defined, i.e., the operator pm,n does not have the 
discrete spectrum. Nevertheless, this continuation provides us a very interesting result. In 
particular, we can see the appearance of the non-perturbative poles analytically. 


5.2 Analytic continuation 

Looking at the differential operators (B.2), one notices that they drastically simplify if 
setting n = —1 (or m = —1) or n = —m. Let us first consider the case of n = —1. The 
computation for n = —m is almost same. In the case of n = —1, the differential operators 
in (B.2) reduce to 


V, 


( 1 ) 


' = —<9^ 

-1 24 


V, 


( 2 ) 

m,— 1 


^9^ 

5760 


V. 


(3) 

m,— 1 


—a® 

967680 ^ 


(5.14) 


The result does not depend on m. Quite interestingly, the same operators also appear in the 
quantum mechanical system with inverted harmonic potential in the computation of the 
quantum periods in the c = I string [49] . Currently, we do not have a clear understanding 
of this agreement because the setups look quite different. The WKB expansion (5.10) also 
simplifies: 


-7 WKB 


(«) 


h 


1 H-H- 

24 5760 


r + 


31s® 

967680 


+ 0{h^) 


(5.15) 


where the classical part is given by 


7(0) 


{s) 


-CSC 

2s 



(5.16) 


It is an easy guess work to find the all-order result to reproduce (5.15). The result is the 
following 


WKB/ X 
m,—l W/ 




h 


hs 

yCSC 



1 /7rs\ 

— - CSC - CSC 

4 \m J 



(5.17) 


We have to note that this zeta function has poles in Res > 0 unlike the ABJM spectral 
zeta. This is of course because of the analytic continuation to n < 0. As mentioned 
previously, for n < 0, the spectral problem is no longer well-defined. Thus we cannot 
apply the argument in the previous section to show that the spectral zeta has no poles in 
Re s > 0. However, as we will see below, this zeta function gives a remarkable result. In 


- 20 - 











the following, we forget a physical interpretation in /r —>■ — oo, and concentrate on the large 
// behavior. It is easy to see that (5.17) has non-perturbative poles at s = 27r^/h {i G Z). 
We stress that these poles appear after the all-order resummation of the WKB expansion. 
We formally plug the result (5.17) into the Mellin-Barnes integral (1.7): 

ds TT , , /TTS' 

■ csc( 7 rs) CSC 


r 

— 1 (/^) ^) — J 


27ri 4s 


m 


CSC 


Hs 




(5.18) 


Now it is clear that the integrand has three-type of poles: two kinds of perturbative 
poles at s = 0, ±1, ±2,... and at s = ±m, ±2m,... and non-perturbative poles at s = 
zb27r/h, ±47r/^,.... In the large // limit, we obtain 

'7r{m + m~^) mh' 




m 




where and are the contributions from the two kinds of the perturbative poles: 


127r;i 
M2,II 


i2h 


+ 


487r 


(5.19) 


OO 


1=1 


= -Y. j; (^) (f) 


TT 

' s=-i As 

e=i 


hs 


^SfJ. 


7 

Ur 


M2,II 
m, — l 


\m J \ 2 
— Res -— csc(7rs) esc 

^ s=-lm As 

e=i 

^ (“1)^ r e \ (-Iran 

^ csc(7r£m) CSC ( j e 

e=i ^ ^ 


/TTS 

V m 


/ hs 


(5.20) 


^Sfl 


Also, is the contribution from the non-perturbative poles: 


WS \—^ T-^ / \ ( I Its 

Jm _l = — / Kes — CSC TTS CSC — esc — 
’ ^ s=- 2 -kIIYi 4s \m) V 2 


( hs 




“ (- 1 )^ 

^csc 


i=i 


Al 


27r7 


CSC 


2^ 

mh 


(5.21) 


e 


These are the complete large /i expansions in the case of n = —1. 

Next, let us proceed to the case of n = —m. In this case, the differential operators 
reduce to 


p(') = —92 

m.—m r> A 


- . 

m-m 24 m-m 5700 ^ m-m 907080 ^ 


V. 


(3) _ Slm*^ 


dl. 


The conjectural all-order spectral zeta function is thus given by 

' mhs' 


^WKB / \ 1 \ ’ 

= --csc(7rms)csc ( - 


Repeating the same computation above, one finally gets 

^7r(l -|- m“2) h 




h 

12TTm?h 


I2h 


I _J_ I .. _|_ /7'M2,I I /^M2,II I /T-WS 

' 48 yr ' ^ ' 'Um,—m ' 'Um,—m ' 


(5.22) 


(5.23) 


(5.24) 
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where 


q-M2,l 
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m 
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'^m,—m 
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£=1 


oo 
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OO 
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(- 1 )^ 

4£ 

(-1)^ 

u 

(- 1 ) ^ 
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csc{7rim) CSC 


imh\ _ 
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■ CSC 
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CSC ( —^ CSC ( ^ ) e ’T.^, 


CSC 


2TT^e 

mh 


0 mh ^, 


(5.25) 


Resolved conifold limit. To make contact with the topological string result in the 
literature, we further take the limit m —>■ 1. In this limit, one should be careful about the 
membrane instanton corrections since both fijygj-gg in m —>■ 1. As seen 

in section 3, such divergences are, however, canceled by each other. The sum of these two 
contributions is finite. After a simple computation, one obtains 


where 


^M2 i T / ^IVlZji , .. / .yiviz.i I .ylViZjli \ 

' m—’ ’ m—' ' 


m —>-1 

OO 


m —>-1 


+ci{h))Q 


1=1 






th /tt , 

ycot -1+1 


There is no limit problem for the worldsheet instanton correction: 


(5.26) 


(5.27) 


^w(/^,^) = E 


1=1 


(-1)' 

it 


CSC 


27r^ 
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e ft 


(5.28) 


Let us return to the geometrical meaning of this limit. It is known that the mirror curve 
(5.1) with (m,n) = (1, —1) describes the resolved conifold (see [59] for example). Therefore 
these results should be compared with the topological string free energy on the resolved 
conifold. In the next subsection, we will see that the results (5.26) and (5.28) are indeed re¬ 
produced by the refined topological string on the resolved conifold, following the procedure 
in [28] . This test is a strong evidence of our proposal that the perturbative resummation of 
Z{s) captures not only the perturbative corrections to J(^) but also the non-perturbative 
corrections in h. 


5.3 Comparison with the topological string on the resolved conifold 

In this subsection, we compare the results (5.26) and (5.28) with the topological string on 
the resolved conifold. We start with the free energy of the refined topological string on the 
resolved conifold with two parameters (€i,e2). 


F{ei,e2;Q) 


OO 
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£^qi/2 _ 


t-^/2) ’ 


(5.29) 
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where 


g = e"i, t = (5.30) 

We first see that the worldsheet instanton correction (5.28) is reproduced from the unrehned 
topological string free energy, i.e., 


ei = -€2 i.e. q = t = qs, Q = Qs- (5.31) 

In this slice, the free energy is given by 

Ftop{qs;Qs) = 112 ^ -1/2.2' (5.32) 

£=i qs ) 

Following [28], we identify the parameters 

qs = e2-\ Qs = e^+-\ (5.33) 

where the string coupling A and the Kahler modulus T are related to h and ^ by 

A = |. T=-A, = -?|ii. (5.34) 

Then the free energy is finally given by 

(~ 1 )^ 2 /27r^.£\ _27r£ 

Ftop{qs;Qs) = J2^l^^^(^~)^ " ■ 

e=i ^ ^ 

This is in perfect agreement with (5.28). 

Next we consider the so-called Nekrasov-Shatashvili (NS) limit [60]: e 2 —)• 0. In this 
limit, the free energy reduces to 


FNs{q]Q) ■■= lim 62 ^( 61 , 62 ; Q) = V 
£2^0 


1 




e=i 


£2 q £/2 _ q - e /2 ' 


(5.36) 


Following [28] again, we identify the parameters as follows: 

q = e~ = Q = e~ = (5.37) 

where A and T is the same parameters above. Then the free energy is written as 


N 00 IT 

I T\ 1 A e^ 
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(5.38) 


Now let us apply the procedure in [28]. This procedure claims that the membrane instanton 
correction to the free energy is constructed from the free energy in the NS limit, 
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Therefore one finally gets 


Fm2 





e 




(5.40) 


This also agrees with (5.26), perfectly. Obviously, these agreements are highly non-trivial. 

Let us explain a significance of these results. In [28], a non-perturbative completion of 
the topological string free energy was proposed. This proposal claims that the sum 


Ftop{qs-,Qs) + FM2{q',Q) (5.41) 

should compute the complete unrefined topological string free energy, including all the 
non-perturbative corrections. On one hand, the standard large N limit of the unrefined 
topological string corresponds to A —>■ 0 or —)• oo with T or n/h held fixed in (5.34). In 
this limit, the genus expansion or its Gopakumar-Vafa resummation captures the first term 
in (5.41). The second term contributes as a non-perturbative correction of order (see 

(5.37)). In particular, the Gopakumar-Vafa resummation (5.35) has an infinite number of 
poles on the real axis of h. These poles are precisely canceled by the membrane instanton 
correction (5.40) as shown for general backgrounds in [28]. As a consequence, the sum (5.41) 
is always well-defined for any h or A. On the other hand, starting with the quantized mirror 
curve (5.3), a natural parameter is of course the Planck constant h. The semi-classical limit 
is —>■ 0. In this picture, the second term in (5.41) is naturally computed. The remarkable 
fact is that this quantum system is just described by the refined topological strings in 
the NS limit [49, 60]. The hrst term appears as a quantum mechanical non-perturbative 
correction [29] . Therefore the unrefined topological strings are closely related to the refined 
topological strings in the NS limit via an “S-duality”. We are seeing the same quantity 
(5.41) from different perspectives.^"^ This is a main consequence of [28, 29]. 

Let us recall that in the computation in the previous subsection we considered only 
the perturbative resummation Z^^^{s). The non-perturbative corrections to Jm,n{fJ-) nat¬ 
urally appear as the non-perturbative poles in the integrand of (1.7) after the resumma¬ 
tion. The fact that the results (5.26) and (5.28) completely reproduces (5.41), proposed 
in [28], strongly implies that the perturbative resummation of Z(s) captures all the non- 
perturbative information on J(/u), precisely. 

Remark. The non-perturbative effects in the topological string on the resolved conifold 
have already been studied in a different approach [62]. The analysis in [62] is based on the 
Borel analysis and closely related to the resurgence theory [12]. See [63], for instance, for a 
pedagogical introduction to the resurgence. The Borel analysis shows that the perturbative 
genus expansion of the resolved conifold free energy is an asymptotic series and non-Borel 

similar but more symmetric situation is also found in vortex-antivortex factorization [61]. In the 
perturbative regime of the vortex partition function, the antivortex appears as a non-perturbative correction 
and vice versa. Here, the roles of vortex and antivortex are played by the worldsheet instantons and the 
membrane instantons. 
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summable. It has a non-perturbative ambiguity^^ of the Borel resummation because of the 
Stokes phenomenon. In [62], the discontinuity of the free energy across the Stokes line was 
computed. The discontinuity is given by 
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where is the string coupling and t is the Kahler modulus. This result is understood as 
a Schwinger effect and thus considered to compute the BPS pair-production rate. Here 
we compare this result with the membrane instanton correction Tm 2 ( 9 ; Q)-^^ To do so, we 
first have to identify the parameters. To match the perturbative free energy (5.35) with 
the one in [62], we have to set 

= 27rA = ^, t = -T-7ri= ^-vri. (5.43) 

h h 

Next, we split the sum over m in (5.42) into two parts: 


Disc Tconifold 


i -T + (2m - l)7ri |(_7^+(2m-l)7ri) 

27rA ^ ^ V ^ ly 

i ( -T - (2m - l)7ri y\ |(_'r_( 2 m-i) 7 ri) 

2vrA ^ ^ V ^ 

i=\ m=l ^ ' 


(5.44) 


After perfoming the sum over m, the first term is evaluated as 


i ^ —T + (2m — l)7ri ^ 


27rA 


EE 

i=\ m=l 

OO . 

E 

e=i 




e 


-T+(2m—l)7ri) 


/m 


CSC 


- 


47r^2 I 2 


m 

ig + — cot 


m 

Y 


+ 1 


(5.45) 




This perfectly coincides with the membrane instanton correction (5.40)! Similarly, the 
second term is the same, but with opposite sign, contribution. Note that in deriving (5.42), 
the refined topological string in the NS limit is of course not used at all. It is remarkable 
that these two fairly independent computations precisely agree. However, naively the sum 
of the two contributions in (5.44) totally vanishes. Currently we do not understand a 
definite reason of this cancellation. It is very important to clarify it more clearly. 


6 Conclusions 

In this paper, we proposed a new perspective on non-perturbative effects in the ABJM 
Fermi-gas system. Our starting point is the spectral problem, and it naturally introduces 

a spirit of the resurgence, this ambiguity must be canceled by the other contributions in the complete 
trans-series expansion. However, the ambiguity (or the discontinuity of the lateral Borel resummation) itself 
has important non-perturbative information. 

^®We thank Marcos Marino for suggesting this interesting comparison. 
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the spectral zeta function. The Mellin-Barnes type representation (1.7) describes the grand 
potential from the small k regime to the large k regime. The spectral zeta function plays 
the essential role in this approach. In particular, it is important to understand its pole 
structure. The consistency condition with the known results requires that the spectral zeta 
function must have the non-perturbative poles, as in (1.9). These poles cause the non- 
perturbative corrections to the grand potential. So far, we do not have a systematic way 
to determine its pole structure, but in this paper, we gave a strong evidence that it indeed 
has the non-perturbative pole at s = —4/A:. Quite interestingly, this non-perturbative 
pole appears after resumming the perturbative corrections to Z(s). Therefore we can 

_ 4/x 

naturally explain the appearance of the worldsheet instanton correction e~~ in the grand 
potential from the perturbative analysis of Z{s). It would be significant to explore the 
pole structure in more detail. It is also important to clarify whether Z{s) itself receives 
the non-perturbative corrections or not. 

We also gave another example on the topological string on the resolved conifold. 
The perturbative resummation of the spectral zeta function perfectly reproduces the non- 
perturbative proposal in [28] (and also [62]). It would be interesting to push this analysis 
in more general cases. 

The Mellin-Barnes representation can be widely applied to many other examples. In 
fact, it also works in the topological string analysis formulated in [30, 36, 37], as seen 
for some special cases in this paper, and in circular quiver Chern-Simons matrix models 
[25, 26, 26, 57] including [24, 53, 54]. In these examples, it is hard to perform the sum 
(1.3) even at the classical level. The representation (1.7) allows us to compute the large 
expansion very systematically. We emphasize that this representation should be valid at 
the quantum level. It would be important to clarify the pole structure of the spectral zeta 
function in those models. 

In the ABJM Fermi-gas, the consistency requires the pole structure (1.9). Naively, we 
expect that Z{s) contains a factor 


OO 


n 

£,m=0 


1 

s + 2£+^' 


( 6 . 1 ) 


A natural function with this factor is the double gamma function r 2 (s|a;i,a; 2 ) with oji = 2 
and u }2 = 4/A:. Also, as in (2.22), the grand potential is simplified by introducing the 
effective chemical potential (2.20). This simplification seems to imply that the poles of 
Z(s) are factorized into 


OO 
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s + 2i 


OO 
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m=l 
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( 6 . 2 ) 


rather than (6.1). It would be nice to consider how such a simplification is understood in 
the Mellin-Barnes representation, and also to investigate a relation between the spectral 
zeta function and the (double) gamma function (or its cousins). We hope to report all of 
these issues in a near future. 
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A Computing the Wigner transform 

In this appendix, we derive (4.5). By definition, the Wigner transform of O is given by 


/ . ipa: / 

Oy^{x^'p)= / dx e n I 2cosh 


x + x'/2\^^^f x — x'/2\^^^ x' 
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The last part is written as 
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Therefore, 


Ow{x,p) = j 


°° dx'dp' Kv-v'W 
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-n' 

— j 2 cosh-, (A.3) 


where we have rescaled the integration variable x' —>■ kx'. Next, we expand the integrand 
around A; = 0. We symbolically expand 




2 X - 2 kx' 
cosh —h smh — 
V 2 4 
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= Cm{x){kx') 
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2m 


(A.4) 


The important fact is that the integral over x' gives the derivative of the delta function: 


dx'_ei(£^(^,)„ ^ (_27ri)”(5W(p_p'). 


(A.5) 


7-00 27r 

Then one can easily perform the integral over p' 

[ —(27ri)^"*(i^^™^(p — p')2 cosh— = (27ridp)^”'2 cosh- = (7ri)^”*2 cosh-. (A. 6 ) 
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Using these results, we finally get 

0\\i{x,p) = 4cosh ^ Cmix){7rik)‘^"' = 4cosh ^ Tcosh^ ^ — sin^ —'j 

m=0 2 


1/2 


(A.7) 


We note that this computation is exact, and 0'w{x,p) contains all the quantum corrections. 
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B Explicit results on differential operators 


B.l ABJM Fermi-gas 

Here we list the forms of differential operators up to n = 4: 


X)(i) 

X)(2) 

X)(3) 

X)(4) 


7r^(l - d^)dl 
96(1+ d^) ’ 

TT\l-d^)dl{lQ + 27d^ + 7dl) 

92160(1 + d^)(3 + d^) 
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+ 12026885,2 + 11358725,^ + 5844805;^ ' 


1 '-7'-7ac\c\ 


+ 31195^ + 1275*). 


(B.l) 


A basic strategy to fix these operators is as follows. We first compute the expansion of 
77^”^ (ft) around k = 0. This can be done by using the formula (4.11) for s = 1,2,... . 
Taking an ansatz of the form of we try to hx unknown parameters to match the first 
several coefficients of If the ansatz is correct, the obtained result must reproduce 

higher coefficients. In this way, one can verify the obtained operator up to any desired 
order. 


B.2 Topological strings for three-term operator 

For the spectral problem considered in section 5, the differential operators are given by 
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967680(m + n + 1) 

As noted in section 5, these operators drastically simplify in m —>■ — 1 or n —>■ — 1 or 
m + n —)■ 0. 
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